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Abstract 

This paper presents an iterated local search for the hxed-charge uncapaci¬ 
tated network design problem with user-optimal flow (FCNDP-UOF), which 
concerns routing multiple commodities from its origin to its destination by 
designing a network through selecting arcs, with an objective of minimizing 
the sum of the hxed costs of the selected arcs plus the sum of variable costs 
associated to the flows on each arc. Besides that, since the FCNDP-UOF is 
a bilevel problem, each commodity has to be transported through a shortest 
path, concerning the edges length, in the built network. The proposed algo¬ 
rithm generate a initial solution using a variable hxing heuristic. Then a local 
branching strategy is applied to improve the quality of the solution. At last, 
an efficient perturbation strategy is presented to perform cycle-based moves 
to explore different parts of the solution space. Computational experiments 
shows that the proposed solution method consistently produces high-quality 
solutions in reasonable computational times. 
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1. INTRODUCTION 

Due to the continuous development of society, increasing quantities of 
commodities have to be transported in large urban centers. Therefore, net¬ 
work design problems arise as tools to support decision-making, aiming to 
meet the need of finding efficient ways to perform the transportation of each 
commodity from its origin to its destination. In the Fixed Charge Network 
Design Problem (FCNDP), a subset of edges is selected from a graph, in 
such a way that a given set of commodities can be transported from their 
origins to their destinations. The main objective is to minimize the sum of 
the fixed costs (due to selected edges) and variable costs (depending on the 
flow of goods on the edges). In addition, fixed and variable costs can be 
represented by linear functions and arcs are not capacitated. Belonging to 
a large class of network design problems, the FCNDP has several variations 
such as shortest path problem, minimum spanning tree problem, vehicle rout¬ 
ing problem, traveling salesman problem and Steiner problem in graph |23] . 
For generic network design problem, such as FCNDP, numerous applications 
can be found mElEH], thus, mathematical formulations for the problem may 
also represent several other problems, like problems of communication, trans¬ 
portation, sewage systems and resource planning. It also appears in other 
contexts, such as flexible production systems [20] and automated manufac¬ 
turing systems [15]. Finally, network design problems arise in many vehicle 
fleet applications that do not involve the construction of physical facilities, 
but rather model decision problems such as sending a vehicle through a road 
or not [23l EH] . 

This work addresses a specific variation of FCNDP, called Fixed-Charge 
Uncapacitated Network Design Problem with User-optimal Flows (FCNDP- 
UOF), which consists of adding multiple shortest path problems to the origi¬ 
nal problem. The FCNDP-UOF involves two distinct agents acting simulta- 
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neously rather than sequentially when making decisions. On the upper level, 
the leader (1^* agent) is in charge of choosing a subset of edges to be opened 
in order to minimize the sum of hxed and variable costs. In response, on the 
lower level, the follower (2"''^ agent) must choose a set of shortest paths in 
the network, through which each commodity will be sent. The effect of an 
agent on the other is indirect: the decision of the follower is affected by the 
network designed on the upper level, while the leader’s decision is affected by 
variable costs imposed by the routes settled in the lower level. The inclusion 
of shortest path problem constraints in a mixed integer linear programming 
is not straightforward. Difficulties arise both in modeling and designing effi¬ 
cient methods. 

The FCNDP-UOF problem appears in the design of a network for hazardous 
materials transportation [31III1II211I9]. Particularly for this kind of problem, 
the government dehnes a selection of road segments to be opened/closed to 
the transportation of hazardous materials assuming that the shipments in 
the resulting network will be done along shortest paths. In hazardous mate¬ 
rials transportation problems, roads selected to compose the network have no 
costs, but the goverment wants to minimize the population exposure in case 
of an incident during a dangerous-goods transportation. This is a particular 
case of the FCNDP-UOF problem where, from a mathematical point of view, 
the hxed costs are equal to zero. 

Several variants of the FCNDP-UOF can be seen on I3i El inmsi El ini [26] 
and have been treated as part of larger problems in some applications on 
ra. The work presented by Bilheimer and Grey |6] formally dehnes the 
FCNDP-UOF. Both Erkut et ah [ 12 ] and Kara et ah [ 12 ] work focus on 
exact methods, presenting a mathematical formulation and several metrics 
for the hazardous materials transportation problem. At Mauttone et ah 
[2n], not only was presented a diherent model, but also a Tabu Search for the 
FCNDP-UOF. Both, Amaldi et ah [3] and Erkut et ah [11] presented heuris¬ 
tic approaches to deal with the hazardous materials transportation problem. 
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At last, Gonzalez et al. [H], presented an extension of the model proposed 
by Kara and Verter [12] and also a GRASP. 

According to [T8l 130] , the simplest versions of network design problems are 
A^P-hard and even the task of hnding feasible solntions (for problems with 
bndget constraint on the hxed cost) is extremely complex [HT]. Therefore, 
henristics methods are presented as a good alternative in the search for good 
solntions. Knowing that, this work proposes an Iterated Local Search Ell 
for the FCNDP-UOF 

This text is organized as follows. In Section 2, we start by describing the 
problem followed by a bi-level and an one-level formnlation, presented on 
|26j . Then in Section 3 we present onr solntion approach. Section 4 reports 
on onr compntational experiments. At last, in Section 5 the conclnsion and 
fntnre works are presented. 

2. GENERAL DESCRIPTION OF FCNDP-UOF 

In this section we describe the problem and present a bi-level and an 
one-level formnlation for the FGNDP-UOF proposed respectively by [91 [26] 
for the FGNDP-UOF. 

The basic strnctnres to create a network are a set of nodes V that represents 
the facilities and a set of nncapacitated and nndirected edges E representing 
the connection between installations. Fnrthermore, the set K is the set of 
commodities to be transported over the network, and these commodities may 
represent physical goods as raw material for indnstry, hazardous material or 
even people. Each commodity k & K, has a flow to be delivered through 
a shortest path between its source o{k) and its destination d{k). The for¬ 
mulation presented here works with variants presenting commodities with 
multiple origins and destinations, and for treating such a case, it is sufficient 
to consider that for each pair {o{k), d{k)), there is a new commodity resulting 
from the dissociation of one into several commodities. 
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2.1. Mathematical Formulation 

This subsection presents a few definitions in order to make easier the un¬ 
derstanding of the problem. 

The model for FCNDP-UOF has two types of variables, one for the construc¬ 
tion of the network and another representing the flow. Let yij be a binary 
variable, we have that i/ij = 1 if the edge [i,j] is chosen as part of the net¬ 
work and i/ij = 0 otherwise. In this case, denotes the commodity Fs 
flow through the arc {i,j). Although the edges have no direction, they may 
be referred to as arcs, because each commodity flow is directed. Treating 
y = hjij) respectively, as vectors of active edge and flow 

variables, mixed integer programming formulations can be elaborated. 

List of Symbols 

V Set of nodes. 

E Set of admissible edges. 

K Set of commodities. 

Set of arcs obtained by bidirecting the edges in E. 

Q Associated graph G(V, E). 

5f Set of all arcs leaving node i. 

6~ Set of all arcs arriving at node i. 

Ca Length of the arc a. 

e(a) Edge e related to the arc a. 
o{k) Origin node for commodity k. 
d{k) Destiny node for commodity k. 
g^j Variable cost of transporting commodity 
k through the arc (i, j) G . 
fij Fixed cost of opening the edge [i,j] G E. 
yij Indicates whether edge [i,j] belongs in the solution. 
xfj Indicates whether commodity k passes through 
the arc (i,j). 
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2.2. Bi-level Formulation 

In FCNDP-UOF, differently from the basic FCNDP, each commodity 
k & K has to be transported through a shortest path between its origin o{k) 
and its destination d{k), forcing the addition of new constraints to the general 
problem. Besides selecting a subset of E whose sum of hxed and variable 
costs is minimal (leading problem), in this variation, we also have to garantee 
the shortest path constraints for each commodity k E K (follower problem). 
The FCNDP-UOF belongs to the class of TV'P-hard problems and can be 
modeled as a bi-level mixed integer programming problem [9], as follows: 


min Y^feVe+Yl H 44 
eeE k^K 

s.t. j/eG{0,1}, yeeE, (1) 

where is a solution of the problem: 


min 


s.t. 


keK a={i,j)^A^ 


E - 

i: 4='’?. 

Vz e F,Vfc G K, 

(2) 

(z,i)G<5+(i) 




< ye, 


= [i,j] G E,\/k G K, 

(3) 

IV 

o 


V(z,j) G A^,Vfc G K. 

(4) 


where: 


6 ? = 


— 1 ii i = d{k), 
1 iii = o{k), 
0 otherwise. 


According to constraints 0-0 , we can notice that the set of constraints 
([^ ensures that the vector of variables y assume only binary values. In (|^, 
we have flow conservation constraints. Constraints ([^ do not allow flow into 
arcs whose corresponding edges are closed. Finally, (|^ imposes the non¬ 
negativity restriction of the vector of variables x^. An interesting remark is 
that solving the follower problem is equivalent to solving \K\ shortest path 
problems independently. 
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2.3. One-level Formulation 

The FCNDP-UOF can be formulated as a one-level integer programming 
problem replacing the objective function and the constraints dehned by (j^- 
(|^ of the follower problem for its optimality conditions [26]. This can be 
done by applying the fundamental theorem of duality and the complemen¬ 
tary slackness theorem [1], as follows: 


where: 


min ^/e?/e+5I H 44 

e^E k^K (ij)gA® 


i: 4- E 4 = '>t 

Vi G F, Vfc G K, 

(5) 




+ a;h < j/g, 

Ve = [i, j] G E,\/k G K, 

(6) 

_k __k \k ^ „ 

~ — ^e{a) — X-a 

Va = (i,j) G A^,k G K, 

(7) 

{ye - = 0, 

Ve = [i,i] G E,\/k G K, 

(8) 

{Ca - 4 + 4 + K(a))x'lj = 0, 

Va = (i, j) G A^, k G K, 

(9) 

A^>0, 

Ve = [i,j\ G E,k G K, 

(10) 

4 e M, 

Vi G F, Vfc G K, 

(11) 

H 

IV 

o 

V(i,j) G A^,Vfc G K, 

(12) 

ye e {0,1}, 

Ve G E. 

(13) 


r —1 if i 

bi = I 1 if i = o{k), 

I 0 otherwise. 


A disadvantage of this new formulation is the loss of linearity of the model. 
To bypass this problem, a Big-M linearization may be used. After it, one can 
write the model as a one-level mixed integer linear programming problem, as 
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follows: 


min 


s.t. 


^ feVe +Y1 H 
eeE keK 

E 4- E 4 = <>?. 

(ij)G(5+(i) {id)&S (i) 


x% + x’i, < ye, 

^e{a) — 

+ MeVe - MeX% - MeX% < Me, 
Me{a)Xij — 7rf + TT^ + Xe(a) — -^e(a) ” Ca, 
^e > 0, 

Trf e K, 

4 G {0,1}, 
ye G {0,1}, 


Vi e V,yk e K, 

(14) 

Ve = [i,j] G E,yk € K 

(15) 

Va = (i, j) e , k € K, 

(16) 

Ve = [i,j] G E,\/k G K, 

(17) 

Va = (i, j) G A^, k G K, 

(18) 

Ve = [i,j] G E,k G K, 

(19) 

Vi G V,yk G K, 

(20) 

y{i,j) G A^,yk G K, 

(21) 

Ve G E. 

(22) 


where: 



— 1 if i = d{k), 
1 if i = o(fc), 
0 otherwise. 


However, optimality conditions for the problem in the lower level are, in 
fact, the optimality conditions of the shortest path problem and they could 
be expressed in a more compact and efficient way if we consider Bellman’s 
optimality conditions for the shortest path problem [T] and using a simple 
lifting process [22]. 


min ^ feVe + X! 

/ . 



eGU feGif (i 

,1)6 A® 



s.t. ^ 4 - 

E 4 = ^^ 

Vi G F,Vfc G K, 

(23) 

(ij)G5+(i) {i 




+ 4 < 


= [i,j] G E,\/k G K, 

(24) 

4 - 4 - ^e(a) 

ye{a){.^e{a) ^a) 

Va = (z, j) G , k G K, 

(25) 

o' 
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Vfc G K, 

(26) 

o"' 

Al 


Vz G \{d(fc)},Vfc G K, 

(27) 

4 e {0,1}, 


G A^,yk G K, 

(28) 

Ve G {0,1}, 


Ve G E. 

(29) 


where: 



— 1 if i = d{k), 
1 if i = o(fc), 
0 otherwise. 


The variables k E K, i ^ V, represent the shortest distance between 
vertex i and vertex d{k). Then we define that will always be equal zero. 
Assuming that constraints (24), (28) and (29) are satisfied, it is easy to see 
that constraints (25) are equivalent to Bellman’s optimality conditions for 
\K\ pairs (o(fc), d{k)). 


3. SOLUTION APPROACH 

This section focuses on presenting the different methods developed in 
this work. First the Partial Decoupling Heuristic is introduced. Secondly a 
procedure to find a lower bound. After that a variable fixing heuristic that 
uses the previously explained methods. At last a Local Branching (used as 
Local Search) and a Ejection Cycle (used as Pertubation) are shown so a 
Iterated Local Search metaheuristics could be done. 
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3.1. Partial Decoupling Heuristic 

The main idea of total decoupling heuristic for the FCNDP-UOF is dis¬ 
sociating the problem of building a network from the shortest path problem. 
This disintegration, as discussed in m, can provide worst results than when 
addressing both problems simultaneously. To work around this situation, 
the method uses what we call partial decoupling, where certain aspects of 
the follower problem are considered when trying to build a solution to the 
leading problem. 

The Partial Decoupling Heuristic iterativily builds a network and then routes 
each commodity so a feasible solution can be built. In order to build the net¬ 
work the cost /g, e G E, k E K is dehned: 



fe + ax + {l-a) X Ce if Ve = 0, 


(30) 


a X of, + (1 — a) X Ce otherwise. 


Doing that we consider whether the edge is open or not, plus a linear combi¬ 
nation of the variable cost and the length of the edge as the hxed cost. The 
a works as a scaling parameter of the importance of the gfj and Ce values. 
In the beginning of the heuristic a prioritizes the variable cost {gfj), while 
in the end it prioritizes the edge length (cg). It is important to pay atten¬ 
tion that gfj = q^/3ij, where represents the amount of commodity k to be 
transported and flij represents the shipping cost through the edge e = (^, j). 
After building the network, another shortest path algorithm, using the edges 
length (ce) as cost, is applied to take every commodity from its origin o{k) 
to its destination d{k) in the built network. 

In order to put the scaling parameter a in good use, the method repeats 
MaxIterDP times and at each iteration using a different value for a. The 
proposed algorithm is a small variation of the original Partial Decoupling 
Heuristic [13] . The procedure is further explained on Algorithm 
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Algorithm 1: Partial Decoupling Heuristic 

1 Input; 7 , A, Q 

2 Data; MinCost cx), a 1, y 0, a; 0; 

3 begin 


4 

K ^ K- 

5 

for numIterDP ml... MaxIterDP do 

6 


while A 7 ^ 0 do 

7 



A CandidateList{K 

8 



k' Random{K)] 

9 



y ir- DijkstraLeader{f^\ k'); 

10 



A ^ K\{k'}; 

11 


for /c G A do 

12 



X DijkstraFollower{c, k)] 

13 


s ^ {y,x)- 

14 


CloseEdge{s)] 

15 


if Cost{s) < MinCost then 

16 



^best ^ ^5 

17 



MinCost ^ Cost{siyest)''> 

18 


Oi MaxIterDP ’ 

19 


A A, x 0, 1 / 0; 

20 

return Sbest 


To solve the shortest path problem, the partial decoupling heuristic applies 
the Dijkstra algorithm. At the \K\ runs, the function DijkstraLeader solves 
the problem of network construction, then, the shortest path problem is 
solved using the DijkstraFollower function, generating a feasible solution. 
The notation s {y,x) means that the solution s is storing the values of 
the variables y and x that were just dehned by DijkstraLeader and Dijk¬ 
straFollower. Since the function DijskstraLeader can open edges that at the 
end do not have flow, we used the function CloseEdge{) set ye = 0 for every 
xfj = 0, V/c G K. The Randomi) function returns a random element from 
the set passed as a parameter. In order to choose the insertion order of the 
I A'I commodities, a candidate list consisting of a subset of commodities not 
yet routed, whose amount is greater than or equal to 7 % times the largest 
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amount (g^) of the commodities not routed is create through the use of the 
function CandidateList{). 


3.2. LBound Method 

LBound Method is a strategy to probably hnd a stronger lower bound to 
the original problem. In order to do that, the method consists in relaxing 
all variables and at each iteration a subset of y variables are turn into bi¬ 


nary variables of the model (23) - (29). The process repeats until |'0.2|T^|'l 
iterations are done or an integer solution has been found. The number of 
iterations was decided after numerical experiments. Details of the method 
could be seen in Algorithm 


Algorithm 2: LBound 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 


Input; K, Q 
Data: nvbin, cont 0 

begin 

Sinf ■(— LinearRelaxion{)] 

E i — 

repeat 

for e E E do 

if Ue > 0.5 then 

!fc e {0,1}; 

E\{eY 

nvbin ^ nvbin + 1; 

Sinf t— SolveRO; 
cont cont + 1; 

until cont > |'0.2|i?|] or OptEound{sinf) = TRUE or 
nvbin > 0.9|i?|; 
return Sinf 


The function LinearRelaxation{) solves the linear relaxation of the problem 
and returns the solution value. The function SolveRQ solves a relaxation the 
problem with a subset of binary variables. Function OptEoundQ verihes if 
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the solution found by the method is integer or not. It is important to remark 
that the condition nvbin > 0.9|ii^| was never reached. 

3.3. Variable Fixing Heuristic 

The Variable Fixing Heuristic (VFH) start using both the Partial Decou¬ 
pling Heuristic and the LBound method. After applying those two methods, 
the VHF uses a relax and £x strategy to try to find a better solution. Based 
on the Relax and Fix Heuristic m, in this third part, we separate the vari¬ 
ables in two distinct sets. Ni is the set of relaxed variables and N 2 is the 
set of binary variables. Initialy Ni contains all variables, while N 2 is empty. 
The main idea is at each iteration move a subset of the flow variables {x^) 
from Ni to N 2 . At the end of each iteration, if a feasible solution for the 
relaxed model was found, the variables y that are both zero and attend to 
the reduced cost criterion for variable hxing, are hxed as zero. The method 
repeats until all have been moved from Ni to N 2 or the duality gap be¬ 
comes lower than one. 

In order to choose the order of variables to become binary, the procedure 
uses a candidate list. To choose a commodity, an element is randomly se¬ 
lected from a candidate list consisting of the commodities whose amount to 
be transported are greater than or equal to 7 % times the largest amount of 
the commodity whose variables are not set as binary. A pseudo-code of the 
method is presented in Algorithm 
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Algorithm 3: VFH 


1 Input; 7 , A, Q 

2 Data; MinCost oo 


3 begin 


4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 
27 


Sbest PartialDecoupling( 7 , K, Q); 

Sinf LBound{K,Q)] 

MinCost Cost{sbest) ; 

K ^ K] 

if OptFound{sinf) 7 ^ TRUE then 

while K and \sbest — Sm/I > 1 do 
k ■«— CandidateList{K 

e {0,1}; 

s SolveR{MinCost)] 

if A feasible solution for the relaxed model was found then 
for e G A do 

if i/e = 0 and RCVF{ye) = TRUE then 

Ve ^ 0 ; 

if Cost{s) < MinCost and Feas{s) = TRUE then 

^best ^ "S ) 

MinCost ^ Cost{sbest) 

else if Cost{s) > Cost{sinf) and Feas{s) = FALSE 

then 

Sinf ^ S, 

else 
I Exit 
K ^ K\{k} 
return Sbest 
else 

I return Sinf 





The function SolveRQ solves a relaxation of the one level formulation (23)- 
(29) with a subset of binary variables, taking into consideration the primal 
bound MinCost. MinCost is dehned as the current best solution cost. The 
RCVF{) function returns TRUE if the Linear Relaxation cost plus the Re¬ 
duced Cost of He is greater than the current VFH solution. The function 
Feas{) returns true if the solution s passed as parameter is a feasible solu¬ 
tion to the original problem and returns false otherwise. 


3.4- Local Branching 

Introduced by Fiscetti and Lodi [13], the Local Branching (LB) technique 
could be used as a way of improving a given feasible solution. The LB 
makes use of a MIP solver to explore the solution subspaces effectively. The 
procedure can be seen as local search, but the neighborhoods are obtained 
through the introduction of linear inequalities in the MIP model, called local 
branching cuts. More specifically, the LB searches for a local optimum by 
restricting the number of variables, from the feasible solution, whose values 
can be changed. 

Formally speaking, consider a feasible solution of the FCNDP-UOP, s = 


{y,x) G P, where P is the polyhedron formed by (23)-(29). The general idea 
would be adding the LB constraint 


(1 “ 2/e) < A, (31) 

eeE|j/e=0 e&E\y^=l 

where A is a given positive integer parameter, indicating the number of 
variables ?/e, e G F, that are allowed to flip from one to zero and vice versa. 
The strategy used here consists on applying the LB constraint only on y 
variables, leaving variables free of LB constraints. 
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3.5. Ejection Cycle 

To understand the principles below the pertubation presented here, it is 
necessary to get to know a few metrics, developed by |2Z] , to evaluate chains 
in a solution. 

Consider a solution dehned by the variables for each arc a G and 
each commodity k & K and ye for each edge e E E. For each open edge e, 
where |/e = 1 and xh > 0 or > 0 for at least one commodity k, the edge 
inefficiency ratio can be dehned as: 

+ xjj) + fe 

h = -n-; Ve = [*, j] 6 E. (32) 

k&K 

The lower the value of fo, more interesting it is to have edge e in the solution. 
The average inefficiency ratio is dehned as: 

^ ^ CUe 

j _ 

e£E 

With these metrics we can dehne a set of inefficient edges as 

Ar = {e \ ye = 1, h > (34) 

As it can be seen above, the set of inefficient edges contains every edge in 
the solution whose inefficiency ratio is greater than the average inefficiency 
ratio. Our aim is to create a movement that remove hows from some of the 
inefficient edges in set Aj. 

After evaluating the edges it is possible to construct inefficient chains from 
a subset of the inefficient edges. First, an edge is randomly chosen from the 
set Ai of inefficient edges to form a component of the inefficient chain. If the 
current partial inefficient chain extends from node i to node j, then an edge 


(33) 
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(a, i) G Aj or (j, h) G Aj is added to the current chain, where nodes a and b 
are not included in the current chain. Whenever an edge is added to a chain, 
it is deleted from Aj. The process of extending the current chain continues 
until no further extension is possible or until the chain is composed by four 
edges. Unless Aj is empty or contains a single arc, the process iterates with 
a random edge chosen to start a new chain. When the process ends, any 
chains containing a single edge is deleted. This is done in order to decrease 
the number of edges affected at each iteration of the method. 

After constructing a set of inefficient chains, we dehne our movement. The 
movement is dehned analyzing each chain in the set of inefficient chains. 
The key aspect of our pertubation is the re-routing of flow from edges of the 
inefficient chain to other edges of the network. First, a list of commodities 
{Kset) that have a positive flow through at least one edge of the randomly 
selected inefficient chain is formed. After that, the opening cost (/e) of each 
edge in the inefficient chain is set as inhnity. After reassigning the costs, 
every commodity in Kset has its route destroyed and reconstructed by the 
Partial Decoupling Heuristic taking into account the new opening costs. If 
a feasible solution is found the method stops, else, another inefficient chain 
is randomly selected and the process restarts. Algorithm describes our 
Ejection Cycle procedure. 
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Algorithm 4: Ejection Cycle 

1 Input; s, 7 , K, Q 

2 begin 

3 

P ^ PInefChain{s)-, 

4 

s •(— 0 ; 

5 

while P 7 ^ 0 and s is not feasible do 

6 


rchain E- RandomfP)-, 

7 


P \ {rchain}; 

8 


Kset ^ SK{s, rchain); 

9 


s E- PartialDecouphng(^, 7, As'£;'r); 

10 


if Cost{s) < Cost{s) then 

11 


s •(— s; 

12 

return s 


In order to clarify Algorithm it is necessary to define a few things. The 
function PInefChainQ returns the set A/ of inefficient chains in a solution 
s. The function SK{) returns the commodities that have a positive flow in 
solution s through at least one arc of the inefficient chain passed as parameter 
and set the hxed costs of the edges in the rchain as inhnity. The function 
PartialDecouplingO reroutes the commodities in Kset- In order to do that 
the DijkstraLeader is applied for all k G Kset and DijkstraFollower for all 
k E K. To account those changes, now the method PartialDecouplingO 
needs to receive a second parameter which is the set o commodities used in 
DijkstraLeader.Besides that a partial solution for all k E K \ Kset is also 
passed as a parameter. 

3.6. Iterated Local Search 

Developed by Lourengo et ah EH. the Iterated Local Search (ILS) is 
a metaheuristic that applies a local search method repeatedly to a set of 
solutions obtained by perturbing previously visited local optimal solutions. 
The ILS presented here uses as its main components, the VFH, the Local 
Branching and the Ejection Cycle presented in the previously subsections. 
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The methods are applied in a straightforward way. First we ran the VFH to 
get a feasible solution and a lower bound. Secondly we try to improve the 
quality of the previously found solution through applying the Local Branch¬ 
ing and the Ejection Cycle. The algorithm is described in Algorithmic 


Algorithm 5: VFHLB 


3 

4 

5 

6 

7 

8 
9 

10 

11 


1 Input; 7, A, K, g 

2 begin 

s,Sinf ^ VFH(^, iF, 7 ); 
s i — LB(s, Ajj 
UpdateBest(s); 

if \cost{sbest) — cost{sinf)\ > 1 then 
while Stop Criterion=false do 
s <(— EjectionCycle(s, l,K,g); 
s i — LB(s, Ajj 
UpdateBest(s); 
return Sbest 


In the VFHLB, the initial solution and the lower bound are generated by 
the VFH method. Then, the function LB performs the Local Branching as 
a Local Search and the EjectionCycle performs a perturbation. 

4. COMPUTATIONAL RESULTS 

In this section we present computational results for the VHFLB presented 
in the previous section. 

The algorithm was coded in Xpress Mosel using FICO Xpress Optimization 
Suite, on an Intel (R) Core TM i3 - 3250 CPU @ 3.5 GHz computer with 
8 GB of RAM. Computing times are reported in seconds. In order to test 
the performance of the presented heuristic, we used networks data obtained 
from Mauttone, Labbe and Figueiredo [26] . 

In order to calibrate the algorithms we use 60% of our data so parameters 
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overfitting could be avoided and the following StopCriterion, 7 and A values 
were tested: StopCriterion = {10 iterations; 50 iterations; 100 iterations}, 
7 = {0.75,0.85,0.90} and A = {[^l, f^]}.After the tests the pa¬ 

rameters were calibrated as: StopCriterion = 10 iterations, 7 = 0.85 and 

The data used are grouped according to the number of nodes in the graph 
(10, 20, 30), followed by the graph density (0.3, 0.5, 0.8) and hnally the 
amount of different commodities to be transported (5, 10, 15, 20, 30, 45). 
We are comparing the VFHLB results with the results of the GRASP pre¬ 
sented by na, which, to the best of our knowledge, is the best heuristic 
aproach to solve the FCNDP-UOF. For the presented tables, we report the 
best solution {Best Sol) and best time {Best Time) reached by each ap¬ 
proach, the average gap {Avg GAP) and the gap {GAP) using the optimal 
solution. We also reported the average values for time {Avg Time) and for 
solutions {Avg Sol). Finally, it is reported standard deviation values for time 
{Dev Time) and solution {Dev Sol). The results in bold represent that the 
optimum has been found. 
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Avg Sol Avg Time Dev Sol Dev Time Best Sol Best Time Avg GAP GAP Avg Sol Avg Time Dev Time Best Sol Best Time Avg GAP GAP 
10 - 0 . 3 - 5-1 3942.00 1.2870 0.0000 0.0329 3942 1.2561 0.0000 0.0000 3942 0.0070 0.0017 3942 0.0060 0.0000 0.0000 

10 - 0 . 3 - 5-2 4552.00 1.3267 0.0000 0.0172 4552 1.3110 0.0000 0.0000 4552 0.0038 0.0004 4552 0.0030 0.0000 0.0000 
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Table 1: Computational results for GRASP and VFHLB approaches for 0.3 density instances 



Avg Sol Avg Time Dev Sol Dev Time Best Sol Best Time Avg GAP GAP Avg Sol Avg Time Dev Time Best Sol Best Time Avg GAP GAP 
10 - 0 . 5 - 5-1 4360.00 1.8240 0.0000 0.0568 4360 1.8058 0.0000 0.0000 4360 0.0086 0.0005 4360 0.0080 0.0000 0.0000 

10 - 0 . 5 - 5-2 1351.00 1.9186 0.0000 0.0632 1351 1.9110 0.0000 0.0000 1351 0.0092 0.0004 1351 0.0090 0.0000 0.0000 
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Table 2: Computational results for GRASP and VFHLB approaches for 0.5 density instances 




Avg Sol Avg Time Dev Sol Dev Time Best Sol Best Time Avg GAP GAP Avg Sol Avg Time Dev Time Best Sol Best Time Avg GAP GAP 
10 - 0 . 8 - 5-1 4033.83 3.0370 108.4895 0.0314 3^6 3.0249 0.1146 0.1014 3619 0.0142 0.0004 3619 0.0140 0.0000 0.0000 

10 - 0 . 8 - 5-2 3535.68 2.9300 60.9944 0.0883 3480 2.9100 0.0160 0.0000 3480 0.1480 0.0016 3480 0.1460 0.0000 0.0000 
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Table 3: Computational results for GRASP and VFHLB approaches for 0.8 density instances 



In Tables were used 135 instances generated by Mautonne, Labbe 

and Figueiredo |2S], whose results were published by them just for 5 in¬ 
stances. For these instances, the computational results suggest the efficiency 
of VFHLB. On average, the time spent by VFHLB was 2.31 times faster than 
the time spent by GRASP, being 2.954 times faster for 0.3 density networks, 
3.167 times for 0.5 density networks and 0.837 times for 0.8 density networks. 
Also, VFHLB found all optimal solutions, while GRASP found only 44 opti¬ 
mal solutions. Besides that, the VFHLB also improved or equaled GRASP 
results for all 135 instances (91 improvements and 44 draws). 

Another important remark is that, in Tables [T] and VFHLB is faster than 
GRASP, both in the mean of Avg Times and in the mean of Best Times. 
Although VFHLB lose to GRASP in the mean of Avg Times and in the mean 
of Best Times on Table On the other hand, GRASP hnds only 26 % of the 
optimal solutions while, as told before, VFHLB hnds all optimal solutions. 
The experiment also showed that, at least for the instances tested, the or¬ 
der of the commodities set by the candidate list in the VFHLB does not 
change the solution obtained at the end of the algoritm, but does affect the 
computational time. 


4.I. Statistical Analysis 

In order to verify whether or not the differences of mean values obtained 
by the evaluated strategies shown in Tables Tp and are statistically signif¬ 
icant, we employed the Wilcoxon-Mann-Whitney test technique [TB]. This 
test could be applied to compare algorithms with some random features and 
identify if the difference of performance between them is due to randomness. 
According to [TH], this statistical test is used when two independent samples 
are compared and whenever it is necessary to have a statistical test to reject 
the null hypothesis, with a signihcance 6 level (i.e., it is possible to reject the 
null hypothesis with the probability of (1 — 0 x 100%)). For the sake of this 
analysis we considered 9 = 0.01. The hypotheses considered in this test are: 
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• Null Hypothesis (HO): there are no signihcant differences between the 
solutions found by VFHLB and the original method; 

• Alternative Hypothesis (HI); there are significant differences (bilateral 
alternative) between the solutions found by VFHLB and the GRASP. 

Table presents the number of better average solutions found by each strat¬ 
egy, for each group of instances separeted by density. The number of cases 
where the Null Hypothesis was rejected is also shown between parentheses. 


Instance Algorithms 


Groups 

GRASP 

VFHLB 

0.3 

0(0) 

30(29) 

0.5 

0(0) 

34(31) 

0.8 

0(0) 

43(33) 


Table 4: Statistical Analysis of GRASP and DPRFLB 

When comparing GRASP with VFHLB, we notice that almost all differences 
of performance (86.91% of the tests) are statistically significant. We can also 
observe that the VFHLB obtained 100% of the best results. These results 
indicate the superiority of the proposed strategy. 

4-2. Complementary Analysis 

Another way to analyze the behavior of algorithms with random compo¬ 
nents is provided by time-to-target plots (TTT-plots) [2]. These plots show 
the cumulative probability of an algorithm reaching a prefixed target solu¬ 
tion in the indicated running time. In TTT-plots experiment, we sorted out 
the execution times required for each algorithm to reach a solution at least 
as good as a predefined target solution. After that, the i-th sorted running 
time, ti, is associated with a probability pi = and the points Zi = (ti;pi) 
are plotted. 

For these experiments we tested 10 of our largest instances with a medium 
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Probabiliiy 


target (1.22 times the cost of the optimal solution). Firstly we analyze the 
instances with 20 nodes, followed by the analyses of instances with 30 nodes. 




(a) TTTPlot: 20-0.8-30-1 (b) TTTPlot: 20-0.8-30-2 




(c) TTTPlot: 20-0.8-30-3 (d) TTTPlot: 20-0.8-30-4 



(e) TTTPlot: 20-0.8-30-5 


Figure 1: TTT Plot - 20 Nodes Instances 
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After analyzing the behavior of the methods for the selected instances of 20 
nodes, through analysis of the TTTPlot £gures[^to[^ we conclude that the 
proposed strategy outperforms the GRASP, since the cumulative probability 
for VFHLB to find the target in less then 40 seconds is 100 %, while for 
GRASP it is 0 %. 




(a) TTTPlot: 30-0.8-30-1 (b) TTTPlot: 30-0.8-30-2 



(c) TTTPlot: 30-0.8-30-3 (d) TTTPlot: 20-0.8-30-4 



(e) TTTPlot: 30-0.8-30-5 
Figure 2: TTT Plot - 30 Nodes Instances 

After analyzing the behavior of the methods for the selected instances of 30 
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nodes, through analysis of the TTTPlot hgures[^to|^ we conclude that the 
proposed strategy outperforms the GRASP, since the cumulative probability 
for VFHLB to hud the target in less then 180 seconds is 100 %, while for 
GRASP it is 0 %. 

CONCLUSIONS 

We proposed a new algorithm for a variant of the hxed-charge uncapac¬ 
itated network design problem where multiple shortest path problems were 
taken into consideration. In the hrst phase of the algorithm, the VFH is used 
to build a initial solution and hnd a lower bound. In a second moment, a 
Local Branching technique and a pertubation. Ejection Cycle, are applied to 
reduce the solution cost. 

The proposed approach was tested on a set of instances grouped by number 
of nodes, graph density and number of commodities to be transported. Our 
results have shown the efficiency of VFHLB in comparison with the GRASP 
presented in [TT], since the proposed algorithm finds the optimal solution 
for all instances and presents a best average time for the majority of the 
instances (125 out 135). 

As future work, we intend to work on exact approaches as Benders’ Decom¬ 
position and Lagrangian Relaxation since both are very effective for similar 
problems, as could be seen in urn- 
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